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Oh' Abstract. To study the evolution of binary star clusters we have imaged 7 systems in the 
Small Magellanic Cloud with SOAR 4-m telescope using B and V niters. The sample contains 
pairs with well-separated components id < 30 pc) as well as systems that apparently merged, as 
' evidenced by their unusual structures. By employing isochrone fittings to their CMDs we have 
determined reddening, age and metallicity and by fitting King models to their radial stellar 
density profile we have estimated core radius. Disturbances of the density profile are interpreted 
as an evidence of interaction. Circunstances as distances between components and their age 
difference are addressed in terms of the timescales involved to access the physical connection of 
the system. In two cases the age difference is above 50Myr, which suggests chance alignment, 
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, capture or sequential star formation. 
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C$ . 1. Introduction 

Differently from the Galaxy, the low-density environment of the Magellanic Clouds 
y—i • favors the gravitational interaction of single clusters, which form binary (or multiple) 
groupings and eventually will merge or disrupt. Pieces of evidence on this subject are: 
(i) Many binary clusters in the Magellanic Clouds have been identified and recognized 
as true physical systems (Bhatia & Hatzidimitriou 1988, Bhatia et al. 1991); (ii) Deep 
and precise colour-magnitude diagrams (CMDs) for isolated clusters in the Large Magel- 
lanic Cloud revealed multi-population sequences (Mackey et al. 2008). If these facts are 
connected their results indicate that we are witnessing different stages of cluster evolu- 
Q\ [ tion in a scenario where cluster groupings and mergers are not uncommon. Evidence of 
cluster binarity and mergers has been reported for both Clouds which are very rich in 
cluster pairs and multiplets (Dieball et al. 2002, Bica et al. 2008, Carvalho et al. 2008). 
We determine astrophysical and structural parameters of SMC binary clusters by fit- 
ting isochrones to the CMDs and King functions to the radial density profiles with the 
■ purpose of providing constraints on the evolution of such systems. 



2. Observations 

The SOAR optical imager (SOI) mounted in a bent-Cassegrain configuration to the 
4.1-m SOAR telescope (Cerro Pachon, Chile) was employed to observe a sample of 7 
merger candidates. Images in Bessel BV filters were obtained on photometric nights 
(Nov/Oct 2007), using the SOI mini-mosaic of two E2V 2x4k CCDs (1 pixel = 15/Ltm) 



119 



120 Santos Jr. ct al. 



Table 1. Cluster sample and derived parameters. 



Cluster 




CV2000 

(h:m:s) 


^2000 

(":':") 


E(B-V) 


ago 
(Myr) 


Z 


(PC) 


A age 
(Myr) 


separation 
(PC) 


age (lit.) 

(Myr) 


NGC220 
NGC222 




0:40:31 
0:40:44 


-73:24:10 
-73:23:00 


0.15 ± 0.03 
0.15 ± 0.03 


70 ± 10 
70 ± 10 


0.004 
0.004 


2.5 ± 1.0 
1.5 ± 0.5 





28 


100 ± 23" 
65 ± 13 fa 

126 c 
70 - 100 d 
100 ± 23° 
70 ± 7 b 

100 c 
70 - 100 d 


NGC241 
NGC 242 




0:43:33 
0:43:38 


-73:26:25 
-73:26:37 


0.05 ± 0.01 
0.05 ± 0.01 


80 ± 10 
80 ± 10 


0.004 
0.004 


2.0 ± 1.0 





9 


79 ± 18 a 

200 c 
79 ± 18 a 
63 c 


B78 
L51 




0:54:45 
0:54:54 


-72:07:46 
-72:06:46 


0.10 ± 0.02 
0.10 ± 0.02 


45 ± 5 
40 ± 5 


0.004 
0.004 


3.3 ± 2.2 
1.5 ± 0.3 


5 


17 




IC1611 




0:59:48 


-72:20:02 


0.07 ± 0.01 


140 ± 30 


0.002 


2.2 ± 1.1 






158 ± 37 a 
100 ± 20 6 

126 c 
100 ± 20 e 


IC1612W: 

IC 1612E= 


=H 86-186 
=IC 1612 


0:59:57 
1:00:01 


-72:22:24 
-72:22:08 


0.07 ± 0.01 
0.07 ± 0.01 


120 ± 30 
60 ± 10 


0.002 
0.002 


2.5 ± 0.5 
2.0 ± 1.0 


- 
60 


- 
8 


180 ± 20 6 

126 c 
50 ± 24 a 
100 ± 50 6 

100 c 


NGC 376 




1:03:53 


-72:49:34 


0.09 ± 0.02 


50 ± 10 


0.002 


3.0 ± 1.3 


- 


- 


32 ± 7 a 
20 ± 2 
16 c 
25 ± 10" 


K50 




1:04:36 


-72:09:38 


0.03 ± 0.01 


50 ± 10 


0.004 


2.5 ± 1.2 






20 ± 5 a 

8 C 


NGC 422 
IC1641 




1:09:25 
1:09:39 


-71:46:00 
-71:46:07 


0.06 ± 0.01 
0.06 ± 0.01 


110 ± 30 
500 ± 30 


0.004 
0.004 


1.4 ± 0.5 
2.9 ± 1.0 


390 


20 





Notes: 

a Pietrzynski & Udalski (1999), b dc Oliveira et al. (2000), c Chiosi et al. (2006), d Mattcucci et al. (2002), 
" Piatti et al. (2007). 



to cover a 5.2x5.2 (') 2 field of view. The CCDs were binned to 2x2 pixels yielding a 
scale of 0.154 "/pixel. The average seeing was ^0.95" in B and ^0.8" in V. Two 
images in each filter were obtained, with single exposure times of 480 s in B and 195 s 
in V. The spatial scale is 1 "=0.3 pc. The CCD frames were reduced with iraf software. 
STARFINDER (Diolaiti et al. 2000) was used to perform PSF photometry. Table 1 presents 
the cluster sample and includes the parameters derived here as well as literature ones. 
The V image of the interesting non-coeval pair IC 1641/NGC422 is shown in Fig. 1. 



3. RDPs and structural parameters 

We built the radial density profile (RDP) for each merger candidate component by 
employing star counts in circular rings around the cluster center. For the sake of unifor- 
mity, all RDPs are based on rings 5" wide. In order to get an optimized RDP, we limit 
the data to a threshold magnitude for which the difference between the cluster central 
density and the adjacent field density is maximum. The cluster center is tweaked in the 
process. All stars fainter than this threshold are then disregarded in the analysis. We cm- 
ployed two-parameter (central density ao, core radius R c ) King profiles (King 1966) to 
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Figure 1. From upper-left to bottom-right: (i) V image extraction (2.6' x 1.7') of pair IC 1641 
(left) / NGC 422 (right); North is up and East to the left; (ii) The RDP in log scale with 
a 2-parameter King profile fitting superimposed for which R c is indicated; 1 a dispersion for 
the fittings are plotted as dotted lines; (iii) Isochrone fittings to the CMDs of IC 1641 and 
NGC 422; for comparison, the best matched isochrone for IC 1641 is also plotted in the CMD of 
NGC 422 and vice- versa; (iv) Limiting radius (Bica et al. 2008) as a function of age; the dots 
represent SMC clusters from Chiosi et al. (2006) and other symbols identify the present sample 
of binary /merger clusters. 



derive structural properties. The background was determined by fitting a constant to the 
outermost 4 rings, which was then subtracted from the overall stellar density before the 
King profile fitting was performed. Fig. 1 shows the RDP for clusters IC 1641/NGC422. 



4. CMDs and astrophysical parameters 

We have obtained astrophysical parameters by fitting Padova isochrones (Girardi et 
al. 2002) to the V x (B - V) CMDs. In order to minimize the contamination of back- 
ground stellar fields, a circular region defined by a limiting radius (Ru m ) covering the 
cluster central area was selected. Rum is defined as the radius where the cluster stellar 
density begins to stand out from the background one. For all clusters in the sample good 
fittings were obtained for (m-M) = 19.0 ±0.1. Reddening as well as age and metallicity 
corresponding to the best matching isochrones are presented in Table 1. Fig. 1 shows 
isochrone fittings to the observed CMDs of IC 1641 and NGC 422. 
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5. Discussion 

In general, our age estimates agree with CMD based age determinations in the litera- 
ture (see Table 1). The ages of L51/B78 and NGC422/IC1641 have been determined using 
isochrone fittings for the first time. Because of a disturbed surface brightness profile and 
isophotal features, Carvalho et al. (2008) classified K50 and NGC376 as mergers and the 
pairs IC1612E/W and NGC241/242 as true binary clusters. Except for NGC422/IC1641 
and IC1612E/W, all cluster pairings in our sample are nearly coeval and have ages close 
to 10 7 yr, consistent with the time scales of dynamical models. NGC422/IC1641 and 
IC1612E/W have ages that differ by 390 Myr and 60Myr respectively, implying that, if 
they are interacting, a capture process may have occurred. Another possibility is that the 
oldest cluster in the pair triggered the star formation that originated the youngest one 
via supernovae explosions and shell expansion. Such mechanism has also been proposed 
as a plausible explanation by Vallenari et al. (1998) for the LMC pair NGC1903/SL357, 
although their age difference is 550 Myr. 

Concerning the time scales involved, the expected effect of the gravitational interaction 
for clusters having a mass of 10 5 M© and half mass radius of 10 pc, is that they will merge 
in 10 7 yr if their separation is ~25pc (Sugimoto & Makino 1989). Models for less massive 
clusters in the LMC tidal field give a few 10 7 yr for the pair lifetime, with separations 
varying from ~6pc for cluster masses 1O 3 M and ~14pc for 10 4 M Q (Bhatia 1990). 
According to these models, after this time scale either the cluster pair merges due to 
mutual gravitational interaction or disrupts by the galactic tidal field and, in a less 
probable event, by passing giant molecular clouds. 

Clusters resulting from pairs that have merged are expected to show structural prop- 
erties which differ from those of single clusters. The plot of Ru m against age indicates 
that K50, NGC376 and IC1611 have a larger Ru m than the mean locus of single clusters, 
which may be an additional evidence for their merger status (Mackey & Gilmore 2003). 
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